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A  new  method  has  been  demonstrated  in  which  a  single  experiment  is  used  to  fully  define 
the  sputtering  angular  dependence  of  a  given  material.  The  method  subjects  a  circular  rod 
of  test  material  to  a  mono-energetic  and  highly  collimated  ion  beam.  The  eroded  profile  is 
then  measured  using  optical  profilometry;  the  full  sputtering  angular  dependence  curve  is 
then  extracted  using  a  numerical  approach. 

I.  Introduction 

Understanding  the  sputtering  characteristics  of  electric  propulsion  thruster  materials  is  important  for  pre¬ 
dicting  lifetime  and  performance.  Sputtered  material  is  also  a  contamination  concern  as  it  can  be  redeposited 
and  degrade  system  performance. 

The  amount  of  sputtered  material  is  a  function  of  the  energy  (velocity)  and  angle  of  incidence  of  the 
impacting  ions.  Although  analytical  sputtering  profiles  exist  for  some  monatomic  metals,1  the  sputtering 
behavior  of  materials  used  in  EP  applications  is  typically  obtained  experimentally  using  weight  loss  and 
quartz  crystal  microbalance  (QCM)  measurements.  In  each  case,  a  flat  plate  is  eroded  by  an  ion  source 
at  a  certain  energy  and  a  certain  angle  of  incidence.  These  traditional  methods  are  characterized  by  the 
time  consuming  task  of  performing  the  multiple  experimental  runs  necessary  to  achieve  sufficient  angular 
resolution. 

The  method  presented  in  this  paper  was  previously  employed  by  Barrie3  et  al.  to  measure  the  sputtering 
rate  of  an  aluminum  rod  subjected  to  bombardment  by  Xenon  particles  from  a  10  cm  Hall  Effect  Thruster 
(HET).  In  this  paper  we  investigate  the  effect  of  the  non-uniform  beam  on  the  erosion  profile  by  eroding 
three  identical  rods  located  at  various  axial  distances  from  the  thruster  as  well  as  differing  angles  from  the 
thrust  axis  as  seen  in  Figure  1  and  Figure  2.  The  increased  distance  from  the  thruster  results  in  a  more 
uniform  flow  profile  at  the  expense  of  a  lower  erosion  rate. 

II.  Model 

Yamamura  gives  the  angular  sputter  yield  as  a  correction  factor  to  the  normal  sputter  yield  with  Equa¬ 
tion  1  below.  The  normal  sputter  yield,  Y(0),  is  calculated  by  Yamamura’s  method  as  summarized  by 
Cheng.4  /  and  S  are  fit  parameters  tabulated  for  a  limited  number  of  source  and  target  combinations  in 
literature.1  0  is  the  incident  angle  relative  to  the  surface  normal. 
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Where  Mi  and  M2  are  the  source  and  target  atomic  masses  respectivly. 
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Table  1  lists  the  required  constants  needed  to  calculate  /  and  E  and  the  angular  dependence  curve  for 
Xenon  sputtering  Al. 


Table  1.  Published  constants  for  Xe  sputtering  Al 


Parameter 

Al 

Mi  amu 

131.29 

M2  amu 

26.982 

Z i,  source  atomic  number 

54 

Z2,  target  atomic  number 

13 

fs 

1.8 

a,  screening  radius 

0.1052 

Ro,  average  lattice  constant 

2.56 

Us,  sublimation  energy  [eV] 

3.39 

Using  the  parameters  above  and  Equations  2-8  above  one  arrives  at  an  /  and  E  of  7.84  and  4.01  respec¬ 
tively  using  the  expected  incident  energy  of  220  eV.  Equation  1  is  plotted  in  Figure  5  with  these  values  and 
a  normal  sputtering  yield  of  unity. 

Since  the  angular  sputter  yield  depends  on  angle  of  incidence  and  incoming  energy,  a  circular  cylinder 
was  placed  in  a  uniform  and  mono-energetic  flow.  When  a  circular  cylinder  is  subjected  to  a  uniform  flow 
perpendicular  to  its  axis,  the  angles  of  incidence  from  90°  to  0°  are  achieved  simultaneously.  One  would 
expect  the  greatest  initial  erosion  to  occur  at  the  angle  of  optimal  erosion,  9opt ,  which  is  in  the  vicinty  of 
60°  off  the  surface  normal.  Furthermore,  the  Yamamura  model  predicts  no  erosion  when  the  flow  is  tangent 
to  the  surface  (9  =  90°)  and  simply  normal  erosion  when  the  flow  is  perpendicular  to  the  surface  (9  =  0°). 
Since  we  cannot  measure  erosion  after  a  very  short  time,  the  surface  is  eroded  over  a  long  period  of  time 
resulting  in  the  surface  normals  changing  significantly.  A  numerical  scheme  was  implemented  to  simulate 
the  changing  surface  normals  and  allowed  for  the  proper  determination  of  /  and  E. 
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Figure  1.  Placement  of  the  rods  in  the  chamber  relative  to  the  thruster,  firing  from  R  =  0  through  a  =  0° 
parallel  to  the  page. 


III.  Results 

In  an  effort  to  subject  the  rods  to  a  more  uniform  and  mono-energetic  flow  of  Xenon  ions,  the  rods  were 
placed  at  a  greater  distance  from  the  source  compared  to  the  previous  experiment  and  at  various  angles  from 
the  thrust  axis.  Locations  beyond  15  degrees  off  the  thust  axis  and  30  cm  from  the  source,  plasma  densities 
drop  off  quickly  reslting  in  lower  erosion  rates.  Figure  1  shows  the  rod  positions  relative  to  the  thruster. 

After  firing  for  135  hours  the  rods  where  removed  from  the  vacuum  chamber,  exposed  to  the  atmosphere 
and  their  surface  profiles  were  measured  using  an  optical  profilometer.  Ten  circumferencial  scans  were 
performed  on  each  rod  across  a  1  mm  swath  centered  at  the  height  of  the  thruster.  Each  of  the  ten  scans 
was  averaged  for  each  rod  and  this  data  was  numerically  smoothed  to  128  points.  Figure  3  shows  the  64 
smoothed  points  from  each  rod  that  experienced  erosion.  The  uneroded  profile  is  plotted  in  the  same  figure 
and  has  a  radius  of  3.98  mm.  Rod  1  obviously  eroded  the  most  because  it  was  the  closest  and  was  placed 
on  the  thrust  axis  whereas  rods  2  and  3  eroded  significantly  less  but  in  a  similar  way  due  to  similar  plasma 
densities. 

The  depth  of  erosion  was  calculated  and  plotted  in  Figure  4  and  wings  of  higher  erosion  are  apparent 
as  sputter  rates  are  higher  at  incident  angles  beyond  normal.  According  to  published  fit  parameters  and 
Yamamura’s  model,  initial  erosion  depth  should  be  the  highest  around  60°,  this  was  not  observed  because 
the  surface  was  eroded  to  a  measureable  depth  and  consequently,  the  surface  normals  changed  significantly. 
As  the  surface  normals  changed,  the  asspumtion  of  a  uniform  flow  striking  a  circular  cylinder  is  not  valid 
and  the  surface  erodes  in  a  complex  way. 

The  fit  parameters  were  extracted  from  the  eroded  profiles  in  Figure  3  using  a  numerical  simulation  which 
applied  Equation  1  on  a  meshed  2D  circular  model  of  the  uneroded  rod.  In  each  simulation,  the  rod  was 
eroded  in  very  small  steps  according  to  Yamamura’s  model  until  the  point  of  normal  sputter  yield  9  =  0° 
reached  the  respective  measured  normal  sputter  depth.  This  effectively  normalizes  the  simulations  in  time 
and  erosion  rates  need  not  be  considered.  An  optimization  routine  was  used  to  determine  /  and  E  values 
which  yielded  best  agreement  with  the  measured  eroded  profiles  and  are  found  in  Table  2. 

The  sputter  yield  angular  dependence  curve  is  plotted  in  Figure  5  for  each  rod  using  the  experimental 
fit  parameters  along  with  with  published  results;  the  plots  were  normalized  using  a  normal  sputter  yield  of 
unity.  It  is  immediatly  apparent  the  experimental  curves  predict  lower  sputter  yields  than  the  published 
curve.  This  could  be  due  to  the  fact  that  the  ion  flow  was  not  truly  mono-energetic  nor  uniform  as  was 
assumed  in  the  experiment.  The  angular  dependece  was  not  captured  with  rod  1  due  as  cross-flow  from  the 
thruster  exit  resulted  in  a  wide  range  of  normal  incidence  upon  the  rod.  Rods  2  and  3  captured  more  of  the 
angular  dependence  due  to  more  ideal  flow  characteristics  at  their  locations  in  plume.  Rods  2  and  3  over 
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Figure  2.  Experimental  setup  and  clearly  visible  erosion 
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Figure  3.  Eroded  profiles  of  the  three  rods  compared  to  the  original  rod. 


and  under  predict  the  angle  of  optimal  erosion  respectively,  this  could  be  due  to  slight  thruster  misalignment 
as  each  rod  was  on  either  side  of  the  thrust  axis.  Figure  6  compares  the  measured  eroded  profiles  with 
the  expected  profile  calculated  numerically  using  published  fit  parameters,  with  •!>  indicating  the  angular 
position  along  the  rod.  In  each  case,  the  simulation  showed  a  characteristic  maximum  erosion  yield  at  an 
angle  of  incidence  of  nearly  60°  as  expected.  Experimenal  results  show  good  agreement  for  shallow  angles 
of  incidence.  As  the  angle  on  incidence  increases;  however,  the  expected  erosion  yield  was  not  observed. 

Figure  7  compares  measured  eroded  profiles  with  their  respective  simulations  using  values  of  /  and  E  from 
the  optimizer  routine.  Figure  8  compares  the  simulated  eroded  profiles  for  each  rod  using  their  respective 
experimental  values  of  /  and  E  with  that  of  the  published  published  values.  In  each  case,  the  surface  mesh 
was  eroded  until  the  normal  sputter  yield  matched  with  experimental  measurements. 

IV.  Conclusion  and  Future  Work 

This  work  was  a  continuation  of  the  experiment  by  Barrie3  in  which  a  new  method  of  determining  the 
sputtering  angular  dependence  characteristics  of  a  material  was  studied.  Current  methods  of  determining 
the  sputtering  characteristics  of  a  material  involve  multiple  time  consuming  steps.  This  paper  further  studies 
the  feasability  of  determining  the  sputter  yield  angular  dependence  in  a  single  measurement. 
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Figure  4.  Erosion  depth  for  the  three  rods  as  a  function  of  angular  position. 


Table  2.  Extracted  /  and  E  for  the  three  rods  and  published  values 


/ 

E 

Rod  1 

0.8 

0.3 

Rod  2 

2.1 

0.9 

Rod  3 

8.9 

5.7 

Published 

7.84 

4.01 

Three  ground  aluminum  circular  cylinders  were  subjected  to  the  plume  of  a  HET  thruster  at  differing 
off  thrust  axis  angles  and  distances  from  the  plume  source.  The  cylinders  were  bombarded  for  135  hours 
and  visible  erosion  occured.  The  eroded  surfaces  were  measured  using  an  optical  profilometer  and  com¬ 
pared  to  their  respective  unerodecl  surfaces.  Optimal  fit  parameters  /  and  E  for  use  in  Yamamura’s  sputter 
yield  angular  dependence  equation  were  then  extracted  from  the  measured  profiles  using  an  erosion  simu¬ 
lator  and  an  optimization  routine.  The  extracted  fit  parameters  yielded  close  agreement  with  experimental 
measurement  s . 

The  fit  parameters  calculated  in  this  experiment  did  not  match  closely  with  published  values,  but  were 
in  better  agreement  than  in  the  first  attempt.  Much  of  the  sputter  yield  angular  dependence  was  captured 
but  not  to  the  correct  scale.  By  moving  the  rods  farther  away  from  the  thruster  exit  and  at  greater  angles 
away  from  the  thrust  axis,  a  more  uniform  and  mono-energetic  flow  was  achieved.  The  more  uniform  flow  in 
this  experiment  subjected  the  rods  to  a  flow  closer  to  that  which  was  assumed.  The  most  likely  reason  for 
the  discrepancy  among  the  fit  parameters  is  the  nature  of  the  HET  plume  structure.  Even  at  the  distances 
the  rods  were  placed,  cross-flow  from  the  thruster  exit  introduced  a  certain  degree  of  flow  non-uniformity 
which  altercl  way  the  rods  were  eroded. 

A  significant  increase  in  accuracy  of  the  experimental  fit  parameters  compared  to  past  work  was  achieved 
simply  by  moving  the  rods  to  a  more  uniform  region  of  the  HET  plume;  one  would  expect  better  results  if 
a  rod  were  placed  even  further.  When  using  a  HET  for  this  experiment,  a  compromise  exists  between  flow 
uniformity  and  the  required  time  to  eroded  to  a  measureable  depth.  A  better  source  could  be  used  such  as 
an  ion  thruster  which  would  yield  the  desired  flow  characteristics  this  experiment  requires. 
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Figure  5.  Sputtering  yield  as  a  function  of  incident  angle. 
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Figure  6.  Comparison  between  measured  rod  profiles  and  simulations  using  published  fit  parameter  values 
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Figure  7.  Comparison  between  measured  rod  profiles  and  simulations  using  extracted  /  and  E  values 
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Figure  8.  Comparison  between  simulations  using  published  and  extracted  /  and  E  values 
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